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Abstract. We study a spatially flat Friedmann model containing a pressureless 
perfect fluid (dust) and a scalar field with an unbounded from below potential of the 
form V(ip) = Wo — Vo sinh (^y^Sj^mp^j , where the parameters Wo and Vq are arbitrary 

and k — ^/WttG~n — Mp 1 . The model is integrable and all exact solutions describe 
the recollapsing universe. The behavior of the model near both initial and final points 
of evolution is analyzed. The model is consistent with the observational parameters. 
We single out the exact solution with the present-day values of acceleration parameter 
qo = 0.5 and dark matter density parameter Q p o = 0.3 describing the evolution within 
the time approximately equal to 2iJ ~ 1 . 
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1. Introduction 

Scalar fields play an essential role in modern cosmology. They are attributed to 
inflation models of the early universe and the models describing the present stage of 
the accelerated expansion as well. There is no unique candidate for the potential of the 
minimally coupled scalar field. Typically a potential is a sum of exponentials. Such 
potentials appear quite generically in a large class of theories (multidimensional [I] - 
[Sj, Kaluza-Klein models , supergravity and string/M - theories, see, for instance, 
[7j and references therein) and from reconstruction and other schemes jH] . 

Single exponential potential was extensively studied within the Friedmann model 
using both qualitative methods QJ and exact solutions [TUj. As to the multiple 
exponential potential, it's not studied well yet. Here we consider some multiple potential 
unbounded from below. According to ^Tj a Friedmann universe involving a scalar 
field with an unbounded from below potential exhibits the interesting features such as 
recollapsing even if it's flat. 

In this paper we study the spatially flat Friedmann model containing a pressureless 
perfect fluid (dust) and a scalar field with an unbounded from below simplest multiple 



exponential potential of the form V((p) = Wo — Vq sinh ( Xy/SnGN'p) , where the 



parameters Wo and Vq are arbitrary and A = y 3/2. For this value of A the model 
is integrable and all exact solutions describe the recollapse of the universe within the 
finite time. The time as well as the intermediate behavior of the model crucially depend 
on parameters Wo and Vq. The paper is organized as follows. In section 2 we describe 
the model and obtain the equations of motion. The equations are integrated in section 
3 using methods developed in our previous papers PQ-[3| and the general properties of 
these exact solutions are studied. In section 4 we fit the model to the observational 
parameters. 

2. The model and equations of motion 

We start with the spatially flat Friedmann universe filled by both a self-interacting 
scalar field ip with a potential V(ip) and a separately conserved pressureless perfect fluid 
(dust). The corresponding Einstein equations read 



where H = a/a is the Hubble parameter, p denotes the energy density of dust, overdot 
means derivative with respect to the cosmic time t, k = \/8ttGn = Moreover, we 

have the equation of motion of the scalar field 






(1) 



(2) 



3Hip - V'(cp) 



(3) 
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and the conservation equation for the perfect fluid 

p = -3Hp. (4) 
The last equation gives immediately 



P = Po{~) , (5) 

where the zero subindex means the present time, as usually. In what follows we consider 
the scalar field potential of the form 

V(cp) = Wo - V sinh ( ^3/2^ , (6) 

where Wo and Vo are arbitrary constants. As the system is symmetrical under the 
transformation (p — > —ip, Vq — > — Vo, without a loss of generality we consider only the 
case Vq > 0. 

Now we introduce new variables x and y by the following transformation 



a 
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xy, K(p = a/2/3 log(y/x), x > 0, y > 0. (7) 



Then the set of Eqs.©,© results in 

x = (tul — col)x — 2u>iu>2y, (8) 
y = 2u) x u) 2 x + {u\ - wl)y, (9) 

where we introduced the positive parameters u)\ and u 2 by 

ujI -uj 2 2 = 3/4k 2 W , Zuj^ = 3/4kV . 

Then Eq.fj^), where the presence of p is cancelled by Eq.fl^J), takes the form of the 
following constraint 

xy - {u\ - ujj)xy - u) x u) 2 {x 2 - y 2 ) = 3/AK 2 p a 3 Q . (10) 

The set of Eqs.tJH]),© may be presented in the following complex form 

Z = L0 2 Z, (11) 

where we introduced the complex variable 

z = x + xy 
and the complex parameter 

uo = iO\ + ioj?. 

It is easy to see that Eq. pi|) implies the following complex integral of motion 

z 2 — uj 2 z 2 = const, 

The presence of this complex integral is equivalent to the existence of 2 real integrals 
of motion. One of them Q(z 2 — uj 2 z 2 ) = — 3/2/€ 2 p flo represents the constraint given by 
Eq. f)l()jl and the other integral 9ft(i 2 — uj 2 z 2 ) has an arbitrary value. 
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3. Exact solutions and behavior of the model 



The equations of motion in form Eq. ljllj) are easily integrable. The result is 

_ (~<„ia ( u)i{t-tn)+iu)2(t-t a -5) _|_ -uJi(t-to)-iu)2(t-to-S)\ 



(12) 



where C > 0. The constants to and 5 are arbitrary, C and a obey the relation 



following from Eq.fjlOj). Then substituting x = 3lz and y = into the relations given 
by Eq.(jZJ) one easily gets the explicit expressions for the scale factor a and the scalar 
field <p. 

Let us now consider the general properties of the exact solutions by analyzing the 
corresponding to Eq.()12|) trajectories (orbits) of the moving point on the Cartesian xy 
plane. It follows from the definition of the variables x and y given by Eq.(|7|l that the 
physical segments of a trajectory belong to the angular domain with x > and y > 0. 
Hereafter we show that each segments is of a finite length and its end-points are attached 
either to one coordinate axis or the both axes. From the physical viewpoint it means 
that all solutions describe the universe evolution within a finite time interval. Moreover, 
as the equations of motion are invariant with respect to the time reflection t — *■ — t each 
segment of the trajectory may be passed in both directions. We notice that the constant 
factor exp(«a) leads to the rotation of the trajectory about the origin (0, 0) through the 
angle a. Further we consider only the case when a = taking into account that all 
remaining trajectories may be obtained by rotation. 

If the time variable t is positive and large enough the last term in Eq. (jl2)l is 
negligible {u\ > 0). Then the motion is confined to the repelling logarithmic spiral 
with x = C exp[u>i(t — to)] cos[a;2(t — T )] and y = C exp[uJi(t — to)} sin[u; 2 (t — To)], where 
T = t + 5. As the time t grows the point clockwise rotates around the origin (0, 0) and 
its distance from the origin increases as exp(c<j 1 t). It means that the evolution starts at 
some moment t = To when the spiral intersects the axis x and finishes when it further 
intersects the axis y. The finite time of evolution is about the value T = T[/{2uj2)- In 
this case the zeros of the functions x(t) and y(t) are different. Then near the initial 
point the expansion of the universe can be approximately described by the power-law 
equation a(t) ~ (t — To) 1 / 3 . The similar equation a(t) ~ [T — (t — To)] 1 / 3 approximately 
describes the collapsing near the final point of evolution. 

If t < and \t\ is large enough the first term in Eq. (jl2J) vanishes and the last term 
dominates. In this case the behavior of the model is similar with the only difference: 
the corresponding spiral is attracting. It describes the counterclockwise rotation with 
approaching to the origin of coordinates x and y. 

The form of trajectories for the intermediate values of the time variable t is more 
complicated. To analyze the behavior one needs to superpose both attracting and 
repelling spirals. The result depends on the constant 5. Evidently it's enough to 
consider uj 2 5 G [0, 7r). If 5 = we get x = C cosh[ui(t — t )} cos[uj2(t — t )] an d 
y = Csinh[u;i(t — to)] sm[uj2(t — to)}. The corresponding trajectory with C — 1 is 




Friedmann universe with dust and scalar field 



5 




Figure lc Figure Id 



Figure 1. The corresponding to the exact solutions trajectories on the 
xy plane. All remaining trajectories may be obtain by rotation about the 
origin (0,0). Each solution is defined on some finite time interval which 
corresponds to a lying in the angular domain with x > and y > segment 
of the trajectory. End-points of such segment correspond to big bang and 
big crunch. 

presented on Fig. la. The moving point passes the curve twice in both directions. It's 
not difficult to prove that for 5^0 any trajectory has no selfintersection points. On 
Figs, lb and lc the constant UJ28 ranges from to tt/2 and from 7r/2 to n, respectively. 
For U2S = tt/2 the trajectory passes through the origin center (see Fig. Id). 

This analysis shows that all exact solutions describe the evolution from big bang 
to big crunch within some finite time. Typically this time is about T = 71/(2002), 
however it may be arbitrarily shorter by choosing trajectories obtained after rotation. 
Each solution is defined on some finite time interval which corresponds to a lying in 
the angular domain with x > and y > segment of the trajectory. End-points of 



Friedmann universe with dust and scalar field 



6 



such segment correspond to big bang and big crunch. If end-points of the segment are 
attached to different coordinate axes then the scalar field goes from — oo to +00 or vice 
versa during the evolution. If both end-points are attached to the coordinate axis y, the 
scalar field diverges to +00 as the time tends to the initial or final value. For almost all 
solutions the scale factor is proportional to the time in power 1/3 near the initial and 
the final points of evolution. But if the end-point of the segment is in the origin (0, 0), 
the roots of the functions x(t) and y(t) coincide and the scale factor is proportional to 
the time in power 2/3 near the corresponding point. 



4. Fitting the model to observational parameters 



After this qualitative description of the obtained exact solution (putting the present-day 
value of time t equal to zero) we present the solution to the equations of motion Eq. (|TT|) 
in the form 



(xo + iyo) cosh cut H — ° sinh cut 
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(13) 



with the values x = x(0), y = y(0), x = x(0) and y Q = y(0) expressed in the terms of 
the observational parameters: the Hubble constant H , the present-day values q and 
Qpo of the acceleration parameter q = a'd/a 2 and the density parameter fl p = n 2 p/ (3H 2 ). 
One easily gets these expressions from Eqs.Q,(j2J) and the definition of the variables x 
and y given by Eq.(JJJ) 



£0 



x 



yo 




(14) 
(15) 

(16) 
(17) 



where we denoted 



£0 



9HI 

8LU1LU2 



i(2 + g )-~fy,o 



The upper sign in Eqs. (116)1 . (|I7}) corresponds to the positive value 0o = ^(0), the lower 
sign appears when ip < 0. We notice that the model implies the value (1 — q )/3 — 
to be nonnegative. It can be expressed via the present-day value w^q of the scalar field 
effective barotropic parameter 

p v 2 /2 - V(<p) 

One easily obtain (1 — g )/3 — fl p o/2 = (1 +w (p o)(l — Q p0 )/2 > 0. Therefore the model is 
consistent with the observational data if the present-day value of the scalar field effective 
barotropic parameter is not less than — 1. 



Friedmann universe with dust and scalar field 



7 




Figure 2. We present the exact solution for the following observational 
parameters: qo — 0.5, Q p q = 0.3 and the following model parameters: 
7tFo/(2w 2 ) = 2, cp = 0, <a, > 0. 



The solution presented by Eqs. (fiT?|) - (fT7]) with the scale factor a and the scalar field 
ip obtained from Eq.(J7|) exists in some finite time interval (t 01 , t 02 ) with negative t i an d 
positive to2 where x > and y > 0. Obviously, for given observational parameters H , q 
and fl p0 the values t i an d t 2 and, consequently, the full time (t 02 — t i) of the universe 
evolution depends on parameters u>i and u)2 determining the potential of the model. As 
we already mentioned in section 3, the typical time of evolution may be approximately 
estimated by the value tt/ (2o; 2 ). Then, instead of Ui and u 2 we use further the following 
dimensionless model parameters: the typical time kHo/{2u)2) of the evolution in units of 
Hq 1 ps 14 billions years and the present-day value Kip(0) = K(p = ^2/3\og{\J 1 + Eq— e ) 
of the scalar field in units of M p = k~ x . Evidently, the model parameter uoi may be 
found for given u 2 and (fo from Eq. (jl8|) . On Fig. 2 we present the exact solution for the 
following observational parameters: go = 0.5, Q p o = 0.3, <^o > and the following model 
parameters: nH / (2ui2) — 2, <po = 0. The interval of definition of the solution is turned 
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Figure 3. We introduce the following moments of time: toi ~ big bang, t$i 
- big crunch, £03 ~~ the beginning of the acceleration, io4 ~~ the end of the 
accelerated stage of evolution, i 5 _ the beginning of the recollapsing. 



to be with the following end-points toi ~ —0.6577Hq an d ^02 ~ 1.2222/7(7 • The scale 
factor presented on Fig. 2a may be described in the main order by a ~ (t — ioi) 1 ^ 3 
and a ~ (t — t^) 1 ^ near the initial and final points. The scalar field presented on Fig. 
2b diverges to +00 as t — > £ i or t — > to2- The scalar field potential given by Eq.(JHJ) 
diverges to —00 as t — > t i or t — > t 02 (see Fig. 2c). Fig. 2d shows the domination of 
the scalar field near the initial and final singularities. Comparing Figs. 2c, 2d one easily 
concludes that the kinetic term 2 /2 dominates the potential V(ip) near these points. 
Then w v — > 1 as t —>■ t i or t —>■ t 02 , i.e. the scalar field is like the stiff matter. 

Now we analyze the behaviour of the model for various present-day values <£>o of 
the scalar field. Besides, t i an d t 02 , we introduce the following important moments of 
time: 

to3 - the beginning of the acceleration, 

to4 -the end of the accelerated stage of evolution, 

tos - the beginning of the recollapsing (see Fig. 3). 
Using the exact solution given by Eqs. (13)-(17) with the scale factor and the scalar 
field obtained from Eq. (7) for the parameter q = 0.5, Q p o = 0.3 (</?o > 0) and various 
(po, we numerically calculate the mentioned time values. The result given in Table 1 
shows the following property: the interval of the accelerated expansion as well as the 
full interval of evolution is shorter if the present-day value of the scalar field is greater. 
The corresponding behavior of the scalar factor is presented on Fig. 4. 
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1 


-0,3607 


-0,0834 


0,0539 


0,2847 


0,4516 


1,5 


-0,2184 


-0,0294 


0,0171 


0,1374 


0,2503 


2 


-0,1251 


-0,0073 


0,0049 


0,0638 


0,1349 


2,5 


-0,0687 


-0,0024 


not found 


0,0269 


0,0711 


3 


-0,0368 


not found 


not found 


0,0122 


0,0368 



Table 1. The dependence of the time values tax, to2j io3j ^04 and ios 
explained on Fig. 3 on the present-day value tpo of the scalar field. 




Figure 4. The scale factor for various present-day values ipo of the scalar 
field. 
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5. Conclusions 

We studied the spatially flat Friedmann model containing a pressureless perfect fluid 
(dust) and a minimally coupled scalar field with an unbounded from below potential 



of the form V((p) = Wq — Vbsinh ( X^/StiGnp) , where the parameters Wq and Vq are 



arbitrary. The parameter A is chosen to be J 3/2 in order to integrate the model in the 
explicit form. All exact solutions describe the recollapsing universe. The behavior of the 
model near both initial and final points of evolution is analyzed. Near the singularity 
the scale factor a ~ t 1//3 for almost all solutions. In this case the constraint Eq.Q is 
dominated by the kinetic energy of the scalar field. Such solutions are called kinetic 
dominated 0. There exists a special solution with a ~ 

£2/3 

near the singularity. It 

appears when the constraint Eq.((TJ) is dominated by the density of dust. It's so called 
fluid dominated solution [Hj. The scalar field typically diverges to ±oo as the time tends 
to the initial or final values. However, there is a special solution where the scalar field 
tends to some constant. 

The evolution between big bang and big crunch crucially depends on the model 
parameters U\ and In particular, the time of evolution is determined by these 
parameters. We show that the model may be consistent with the observational 
parameters. We singled out the exact solution with the present-day values of the 
acceleration parameter go — 0.5 and the dark matter density parameter Q p q = 0.3 
describing the evolution within the time approximately equal to 2i/ ~ 1 . 
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